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APPROXIMATIONS IN THE INFINITE POPULATION MODEL 


OF THE ARPANET SATELLITE SYSTEM 


1: Introduction 


In ASS Note 12, а "МатКоуіап" model for the satellite channel with an 
infinite user population* was studied in detail. The relevant questions 
regarding channel throughput rate and packet delays in relation to traffic 
levels and random retransmission delays** were studied. However, such results 
were based upon a very strong assumption (A.4 in ASS Note 12): The traffic 
arriving to the satellite channel is assumed to be a Poisson process. This 
assumption implies two things: (A.4.1) the number of arrivals to a time slot 
is Poisson distributed and (A.4.2) the arrivals to all time slots are inde- 
pendent. Such a strong assumption has invited some critical comments. In this 


note, for a large k we relax assumption (A.4) to 


(A.4') arrivals to different time slots which lie within k slots 


of each other are independent. 


Using (A.4'), we prove that the traffic arriving to a time slot indeed 
has a Poisson distribution as К + о, We also elaborate upon the approximation 
involved in this mathematical model. Lastly, an improvement is made to the 
results of ASS Note 12: q 8 Prob [successful transmission|new packet] will be 
calculated in a fashion which is more consistent with the "Markovian" assump- 
tion of this model. As a result, we obtain slightly more optimistic results. 
The overall behavior of the model is unchanged. Also, all the equations given 
in ASS Note 12 remain valid. Proofs for propositions 1 and 2 in this note are 


* 


*This assumption gives pessimistic results since in this case, every 
packet will interfere with every other packet. 


**We again assume that the random retransmission delay is uniformly distri- 
buted over k slots. 
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given іп Appendix І. In Appendix II, we give some comments оп Bob Metcalf's 
ASS Note 16 and suggest an alternate solution. 


random retransmission delay 








blocked packets 


new arrivals at rate O 


traffic throughput rate - O 


channel 


Fig. 1. The Model 


2. Distribution of Traffic Arrivals to a Time Slot 


Any dependence among arrivals to different time slots is caused by pre- 
vious collisions at the satellite channel. Given two time slots within k 
slots of each other, if КТ < d,* there will be no direct retransmitted 
packets from one slot to the other. Hence, any dependence arises when collided 
packets from the past are retransmitted to the given slots. Such dependence 


is believed to be small and decreases as k increases. 


Define 


c, 8 number of traffic arrivals to the nth time slot 


pf") 2 ртоь{с = i] iow Dil э 


*T is the slot length and d is the round trip propagation delay. For a 
satellite channel d is large and КТ < 4 is true in almost all practical 
situations. However, in the theoretical limit as К + о ме must assume that 
T* 0 if we are to have kT « d. 
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9,2) à 2, бА 


A lim ,(n) 
Р, 7 mo Pi 


иода) 


Q(z) 


Y 9 Y iP, = average channel traffic/slot 
іші 


vi Ё Prob[number of external new arrivals to a time slot = i] 
со 


V(z) А > RA 
i=0 


с В 2, іу, = average throughput/slot 
ізі 


Without loss of generality, let T= 1 


Proposition 1 
Assuming (A.4') and under stationary conditions, Q(z) is given by 
Р k 


да)- |а - 2 + а - +2) | ve) a) 


If the external new arrivals are Poisson distributed, the traffic arrivals 


are Poisson distributed in the limit as k >œ, and 


La Gapee re | (2) 
Lm os ye" (3) 


k>% 
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Eq. (1) is quite general requiring only assumption (A.4'). The external 
new arrivals can have any distribution. An explicit solution for Q(z) cannot 
be obtained except when К >œ or in some uninteresting cases, one of which 


we give here as an example. 
Example k=1 
We have from Eq. (1) 


2 2 2 
Po + Pi 2% Р, Zo» vox 7 (Pg * Pi * Р, zZ + JG *tYztvz t conn 199 


Equating coefficients for the first two terms, we obtain 








ley 
0 
Р, = P 
1 Vo 0 
У 
М 1 
e кси fe 
1 
which imply that Eq. (4) is inconsistent unless v, = 1 - у. Hence, as 


zd 0 
п + © the stationary limit Q(z) does not exist unless Prob[number of new 
arrivals to a slot > 2] = 0. In that case, 


and from Eq. (4) Q(z) = (Po + Pj) (vg E у12) = у 


trivially true. 


o * у12 = V(z) which is 


Se Approximations in the Model 


Since a mathematical model is only as good as its assumptions, in this 
section we shall briefly go over the model approximations. Consider the 
arrival of all packets (new or retransmitted) to the satellite channel as a 
discrete time stochastic process. To give an exact characterization of the 
process, we Une joint probabilities e = ір а, = i», soie $ ЕЯ = ig) 
for all k, пу, ...› Пу, ii, Ге і, > 0. From these we get both the distri- 
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bution of the number of arrivals to any time slot as well as the time correla- 
tion among arrivals to all time slots. A first-order approximation (as done 

in the early ASS Notes) is to assume Poisson distributed arrivals to each time 
slot and independence among arrivals to all time slots. The Poisson assump- 
tion has been shown by the previous proposition to be good when k is large. 
The independence assumption, however, has been weakened to that of "local inde- 
pendence" among arrivals to k adjacent slots. Time correlation is assumed 

to be "Markovian," i.e., the process has a finite memory about its past. Such 

a model is believed to portray fairly accurately the behavior of the real system 


without getting too involved in mathematical analysis. 


4. An Improvement to Results in ASS Note 12 


The "Markovian" time correlation assumption has the added advantage of 
bringing out the effects of the random retransmission delay (assumed to be 
uniformly distributed over k slots). Recall from ASS Note 12, 


q = Prob[successful transmission |new packet] 


4. Prob[successful transmission|previously blocked] 


The effect of k was built into а,› whereas q was taken to be Poi pii 


Since the event [new packet] gives us added information, the effect of k 
should also be built into q to be consistent. In a similar manner to the 


derivation of 4.» ме now see that 
gu 44 e (5) 
where, referring to Fig. 4 in Appendix I and for Т=1 


Prob[no retransmissions from the (п-4-1) slot to slot п] 


qj = 
«YE VY ey 
and 
PME. НИВА 
Tow vr 
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Proposition 2 
ҮТ 


Іп the limit as k +œ, q converges to е. 


D 


5. Final Remarks 


In this note we have examined the approximations involved in our mathe- 
matical model for the infinite population satellite system. We came to the 
conclusion that the Poisson and "Markovian" assumptions will serve our purpose 
in approximating the behavior of the traffic arrival process to the satellite 
channel. 

An improvement was made to give q a new value to reflect the "Markovian" 
property of the model. In Figs. 3 and 4 we plot q and 4. and the channel 
throughput rate 0 as functions of k and Y. Comparing these with the cor- 
responding plots in ASS Note 12, it can be seen that we now have more optimistic 
results for small values of К and ү. Also, all the equations in ASS Note 12 


remain valid with the new value for q. 
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APPENDIX I 


Proof for Proposition 1 
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"Ev 


new packets 


Fig. 4. Traffic Arrivals to a Time Slot 


Let oe = Prob[i retransmissions from slot (n-d-j) to slot n] 
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=. wall Р = = Ya Е 1 (n-d-j) 

= =a > 2, т # 1-2) Р 
(n-d-j) 

P со 

-"——а-5,Б(-{+#) p(n-d-j) 
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Assuming (A.4') and since the external new arrivals аге assumed to be independ- 


ent of the state of the system, 


же 
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"n 
о 


Q, (2) (1 dca je) уа) (1.1) 


j 


Under stationary conditions as п > о, 


[Qc ]* vez) 
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Ia - z) +Q (1 - it 2) V(z) 


where ба) 8 = ei 4. 4-3 (2) 


Q(z) 
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1=0 1=1 
1-10-29 + o(1) 
Q(z) = | .Yü-2, а 5% | V(z) (1.2) 
m Q(z) = ee vb 


Now assuming the external new arrivals to be Poisson distributed, 


v(x) = e 9072) 


йа T (72) o(1-z) 


koo Q(z) 


e 


where Pi = Prob[exactly 1 arrival to the satellite channel in a time slot] 


Prob[a packet obtains successful transmission] 


throughput rate of the channel. 


Under stationary conditions, the channel throughput rate is equal to the 


channel input rate, 


Hence, 


m Q(z) = = which is the z-transform of a Poisson distribution 


with arrival rate = y per slot. An immediate consequence is, 


lim -Y 
kao © = Ру = ye 


Q.E.D. 


ASS NOTE 17 12 


Proof for Proposition 2 


Assume 4 = as we want to show that Eq. (5) holds in the limit as 


k > œ 
Let T= 1 
k k -[q,/(q,+1-q) ЈҮ - 
lim lim Y сүтүк -Y/k ы" 
нз mis = а [a+ Le | E T 


From Proposition 1 in ASS Note 12, we have 


Wen а, = e` given that а а = е“ 


Hence, 
: m MU. CN Р РУТЬ: PNIS | 
as mig a e ИЫ 7e ее )JN 
ве! 
lim 
"ұл LHS 


Q.E.D. 
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APPENDIX II 


Comments on ASS Note 16 by Bob Metcalfe 


In ASS Note 16, with Q defined to be the steady-state time-average of 
the number of terminals with packets ready and X = probability of transmission 
at a terminal with a ready packet,* the following equations are given (in 


Fortran notation): 


W = probability of throughput in a slot 
уле ри а 29 (11.1) 
U = probability of traffic in a slot Г 
=1-({(1-Х) "4 (11.2) 
Xopt = value of X which optimizes throughput rate for a given Q 


1/Q (II.3) 


We discussed these equations with Metcalfe and pointed out that his use 
of the average value Q, was improper since he was dealing with non-linear 
functions of Q. We compared his results with ours (given below) and agreed 
that as X> 0, both his and ours yield the same results for Eqs. (11.1) and 
(11.2); however, Eq. (11.3) must be modified as we show below. 


If we define the random variables 


instantaneous number of terminals with packets ready 


o 
Ш 


w = instantaneous probability of throughput in a slot 


u = instantaneous probability of traffic in a slot 


and observe that if Q, W and U are the instantaneous conditional values of the 
random variables q, w and u, then the above equations are indeed correct. 


*In this case the packet retransmission delay is geometric rather than 
uniformly distributed. 
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However, in this case only the instantaneous value of the throughput rate in 


the next slot is maximized with Хор = 1/4. It is not clear if the steady- 


state throughput rate is also maximized. Also, system overhead will be sig- 
nificant since q has to be measured and the value X 
slot. 


opt adjusted at each 

It is assumed in ASS Note 16 that the slot duration D is long enough 
so that a sender will know the success of a transmission before the start of 
the next slot. This assumption may be valid for a radio channel, but is 
unrealistic for a satellite channel due to the large round trip propagation 
delay. In the latter case the set of terminals with packets ready may be very 
different in 2 adjacent slots and it is clear that we need more than just the 
first moment of q to obtain the optimal value of X. 

We propose the following alternate solution. The value of X is adjusted 
periodically. Assuming that the probability distribution (actually we need 
only the first 2 moments) of q in the current period can be measured and it 
will approximate the probability distribution of q in the next period (this 
assumes the length of a period to be short), we rewrite Eqs. (11.1) and (11.2) as, 


B 27 1 
ЕМІ = gs X P, (II.4) 
N А 
Elu] =1- $ 0-х)" ^ (II.5) 
150 


where N = total number of terminals 
Prob[q = i] 


"Y 
[LU 


Putting эн] = 0, 


* 


we have 


ы i-2 
2s id- iX(- X) Р. = 0 


Assuming that X is so small that we may neglect higher order terms in the 


expansion of (1 - iX)(1 - 0", we have 
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= ifi « Dii nx) P. = 0 


i=0 z 
4-2Х(42- 9) = 0 where q? = E[q"] 
ОРЕ ПИНК) ЕЛИЙ 
РЕ | (11.6) 


2042-9) 24/6 - 1) 


The above is the optimal value of Х for the static optimization of 
E[w] given that we know the first two moments of q. 

From the above discussion, we see that if we know the exact value of q 
(the number of terminals with ready packets), then. X = 1/4 is optimal. How- 
ever, if we do not know the exact value of q, then the use of 8. 5 ШУО 


X = 1/q is not optimal. Instead 


1 


2a /q - 1) 


Xopt ^ 


under the approximations of our expansion. 


